Orbital ordering in cubic LaMnOs from first principles calculations 



H. Zenia,^'^'0 G. A. Gehring/ and W. M. Temmerman^ 

^Department of Physics and Astronomy, University of Sheffield, Sheffield, S3 7RH, UK 
^ Daresbury Laboratory, Daresbury, Warrington WA4 4^D, UK 

We report on first principles Self-Interaction Corrected LSD (SIC-LSD) calculations of electronic 
structure of LaMnOs in the cubic phase. We found a strong tendency to localisation of the Mn eg 
electron and to orbital ordering. We found the ground state to be orbitally ordered with a staggered 
order of ~ and y^ — orbits in one plane and this order is repeated along the third direction. 
The difference in energy with a solution consisting of the ordering of — and 3i/^ — is, 
however, very small. The system is in the insulating A-type antiferromagnetic ordered state in 
both cases. The presence of orbital ordering means breaking of the cubic symmetry and without 
recourse to distortion. The latter may rather be the result of the orbital ordering but the symmetry 
of this ordering is determined by coupling to the lattice. The strong tendency to localisation of 
the eg electron in LaMnOa accounts for the survival of local distortions above the structural phase 
transition temperature. 



I. INTRODUCTION 

Orbital ordering in LaMnOa has always been associ- 
ated with the Jahn- Teller instability of the system as 
a result of the degeneracy of the orbits. Because of 
the crystal field due to the oxygen octahedron the i2g 
orbits are known, from a local picture, to lie lower in 
energy than the ones. The t^g orbits are occupied 
each by one electron and because of strong intraatomic 
Hund exchange interaction the spins of these electrons 
are aligned parallel forming a S" = 3/2 spin which is 
sometimes treated as a classical spin in model calcula- 
tions. In LaMnOa there is one more electron that occu- 
pies an Cg orbit. Because the two eg orbits are degenerate 
the system is unstable towards a distortion which would 
lift the degeneracy. The amount by which the system is 
distorted is then determined by the competition between 
the gain in the electronic energy and the increase in the 
elastic energy of the lattice due to distortion. LaMnOs 
is found in the distorted phase below 780K. Orbital or- 
dering in LaMnOa has been observed by resonant X-ray 
scattering on the Mn-K edge. It was found that this or- 
dering decreases above the Neel temperature 140K and 
disappears above T=780K concomitant with a structural 
phase transitional. 

In 3d transition metal compounds with orbital degen- 
eracy two scenarios are invoked to explain orbital order- 
ing. On the one hand there is superexchange interaction 
between orbitals on different sites. This involves virtual 
transfer of electrons and strong on-site electron-electron 
interaction. On the other hand cooperative JT distor- 
tions or electron-lattice interaction leads to splitting of 
the degenerate orbits and thus to orbital ordering. Al- 
though the ground state ordering of LaMnOa can be ex- 
plained by both mechanisms it is not easy to say which 
is the dominant contribution. This question may sound 
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of little importance so far as LaMnOa is concerned but 
it is important to know the answer because whichever 
is the dominating mechanism will remain more or less 
active once the system is doped. In which case the two 
mechanisms, i.e., electron-lattice interactions or electron- 
electron correlations may lead to different physics for the 
doped systems. 

There have been attempts, by using model calculation, 
at explaining how orbital ordering can occur if one as- 
sumed an antiferromagnetic spin ordering However, 
as mentioned above, temperatures at which the orbital 
ordering sets in are much higher than the Neel temper- 
ature of A-AF spin ordering. Orbital ordering can not, 
therefore, be attributed to spin ordering. In previous 
LSD 0,1113 and model HF [Ell calculations it is found 
that inclusion of distortions is necessary to recover the 
correct A-AF and insulating character of LaMnOa in the 
ground state. The cubic system was found to be both 
metallic and ferromagnetic in the LSD calculations. By 
using Self-Interaction Correction (SIC) to the LSD we 
can allow for the t^g orbitals to localise and form low ly- 
ing semi-core manifold well below the eg orbits 01 . We 
then can compare total energies for different scenarios 
corresponding to localising a particular eg orbit. By do- 
ing so one breaks the cubic symmetry but this is allowed 
if the resulting ground state is lower in energy. As a 
result of the orbital ordering the system will distort in 
order to reduce the electrostatic energy due to the inter- 
action of the oxygen electronic clouds with the lobes of 
the eg orbit, that is directed to them, on neighbouring 
Mn ions. 

The basis states are written: \x) — x'^ — and 
\z) — -^{2z'^ — x"^ — y'^). A composite state can be writ- 

ten as: \9) = cos||z) -I- sin||a;) p. Then the orbital 
state \z) corresponds to 6* = and the state \x) to 9 = tt. 
The orbital ordering of LaMnOa consists of an antiferro 
ordering of two orbits, viz., \±9) — cos ± sin in 
a plane while the same order is repeated along the third 
direction. Until recently it was assumed that 9 = 27r/3. 
But recent ESR [l^ and neutron diffraction mea- 
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TABLE I: Total energies in mRy per formula unit and magnetic moments in fiB of cubic LaMnOs in the FM, A-AFM and 
G-AFM magnetic orderings with several orbital ordering scenarios. Where one orbit only is specified the orbital ordering is 
ferro and for two orbits e.g. — /Sy^ — the ordering is of a C-type with the ordering vector q = -^(1,1,0). The energies 
are given as differences with respect to the energy of the solution corresponding to the experimentally known structure of the 
distorted LaMnOs. 

Configuration | Isd t2g 3z^ - x'^ - y"^ - /3y^ - x^-z^/y^-z'^ - /3z^ - r'^ 

FM 1403 2L4 O TTT 05 ITE 03 

Energy A-AFM 152.0 30.7 7.2 11.4 0.0 -0.6 4.9 

G-AFM 160.9 45.2 9.4 9^5 7^5 7/7_ 8^8 

FM 2l9 Sm 172 Tm 370 sTl 3.68 (3.72) 

Mn mom. A-AFM 2.81 3.14 3.62 3.69 3.70 3.67 3.67 (3.64) 

G-AFM 3.10 3.41 3.60 3.60 3.60 3.61 3.63 (3.57) 



surements have etimated 9 to be 92° and 106° respec- 
tively. Phenomenological superexchange calculations for 
the ground state ordering have also given 9opt 83°, 
"significantly different from 27r/3" 0- Our current cal- 
culations are however limited to the cases of ferro order 
of 6 =0 and tt and antiferro order of ±n/3 and ±27r/3. 



II. CALCULATION DETAILS 

The calculations are performed in the SIC-LSD 
approximation implemented within the LMTO-ASA 
method JJj . The SIC corrects for the spurious in- 
teraction of an electron with itself inherent to the LDA 
approximation to the exchange correlation functional of 
the DFT. It is known however that this energy is im- 
portant only when the electron is localised whereas it 
vanishes for delocalised electrons. This method is used 
to determine whether it is favourable for an electron to 
localise or to be itinerant. This is done by comparing 
the total energies of the system in the presence of the 
two scenarios. The lattice parameter used in the present 
calculation, oq = 7. 434a. u., is the one which gives the ex- 
perimental volume of the real distorted LaMnOa system. 
We have used a minimal basis set consisting of 6s, 5p, 5d 
and 4/ for La, 4s, 4p and 3d for Mn and 2s, 2p, and 3d 
for O. Mn Ap and O 3d were downfoldcd. For the atomic 
sphere radii we used 4.01, 2.49 and 1.83 a.u for La, Mn 
and O respectively. In order to look at different orienta- 
tions of the two orthogonal eg orbitals we used rotations 
of the local axes on the Mn sites. We checked the accu- 
racy of these rotations by comparing the total energies 
of three configurations: all 3z^ — r^, all 3x^ — and all 
3y2 _ ^2 localised in both FM and G-AFM cases because 
these magnetic orderings preserve the cubic symmetry 
and hence the energy should not be dependent on which 
orbit is localised so long as it is the same one on all the 
Mn sites. The energy differences found in this way were 
always less than ImRy per formula unit. 

The calculations were done for a four-formula unit cell. 
The notations of the orbital ordering scenarios are as 
follows: Isd: LDA calculation with no self-interaction 
correction; t2g'. SIC applied to the t2g orbits only on all 



the Mn sites, and in all the other cases one eg orbit is 
localised on top of the t2g ones. The remaining scenarios 
correspond to localising either the same or different orbits 
in the ab plane while preserving the same ordering on the 
second plane along c. Thus we have either ferro or C-type 
antiferro orbital ordering. 



III. RESULTS AND DISCUSSION 

From the total energies of Table we see that the 
ground state corresponds to an orbitally ordered solu- 
tion forming a C-type antiferro-orbital arrangement of 
the x'^ — and — in the ah plane with the same 
ordering repeated along the c— axis. The corresponding 
magnetic ordering is of A-type AFM as found in the dis- 
torted system. This solution is however almost degener- 
ate with the solution with an ordering of 3%^ — and 
3^2 _ ^1 ^ rpj^g energy difference between the two solu- 
tions, 0.6mRy/f.u., is within the accuracy of the calcu- 
lation method (LMTO-ASA). It is then most hkely that 
the true ground state of the cubic system is made up of 
a combination of both solutions. Interactions with the 
neighbouring oxygens are certainly different for the two 
orderings and relaxation of the oxygen positions in the 
real system may favour one of the solutions or a linear 
combination of them. 

We have considered three types of spin order. Ferro- 
magetism and A type antiferromagnetism where the spins 
are parallel in the x-y planes and the planes are stacked 
antiparallel up the z axis and G type antiferromagnetism 
where each spin is antiparallel to all its neighbours. The 
difference in energy between the FM and A-AFM mag- 
netic orderings in the two cases is also very small which 
is consistent with the fact that inter-plane AF exchange 
is much smaller than in-plane FM exchange in agree- 
ment with experiments. Experimental exchange inte- 
grals are obtained from fitting neutron scattering results 
(spin wave dispersion) to a simple Heisenberg Hamilto- 
nian with two exchange integrals acting between nearest 
neighbours. We calculated the exchange constants us- 
ing the convention of Ref Hi]: Ep = (-4Ji - 272)5^, 
Ea-af = (-4Ji + 2J2)^2 Eg-af = (4Ji + 2J2)S^ 
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TABLE II: Magnetic exchange constants in meV obtained from the total energies in Table|I] Ji and J2 are Heisenberg in-plane 
and inter-plane exchange integrals respectively. 



00 scenario 


2 2 
3z — r 


2 2 
X — y 


3a;2 - rVSy^ - 


2 2/2 2 

.T — z /y — jz 


3x'^ - rVSz^ - 




14.96 


-12.93 


51.02 


56.46 


26.53 


4J2S 


-6.12 


-2.04 


-3.40 


-0.68 


-9.52 



for the energies of the FM, A-AFM and G-AFM respec- 
tively. We assumed the value of = 2 for the mag- 
netic moment on Mn ions for all the orderings. The re- 
sults are given, in Tabled for different orbital ordering 
(00) scenarios of the eg orbits as given in Tabled Ex- 
perimentally the two exchange integrals are found to be 
8J1S = 13.36 ± O.lSmeV and 4J2S = -4.84 ± 0.22meV 
for the in-plane and inter-plane coupling respectively 
[1^ HH l . We see then that our calculation overestimates 
the tendency to in-plane ferromagnetism whereas the in- 
terplane exchange is marginally underestimated. How- 
ever it was found in LDA calculations 0] that the first 
neighbour exchange integrals depend dramatically on lat- 
tice distortions. This might explain why our exchange 
constants calculated for the cubic lattice are quantita- 
tively different from the experimental ones which were 
determined for the distorted lattice. Our results are how- 
ever in disagreement with recent model calculations of 
Sikora and Oles who found that for an ordering of 
"0 = 27r/3, often assumed for LaMnOa," the exchange 
constants "are never close to the experiment" . Their cal- 
culated constants are both ferromagnetic which contra- 
dicts the experimental fact that LaMnOa is an A-type 
antiferromagnet. Hence their argument that should in 
fact be different from the assumed 27r/3. 

The widely used Goodenough-Kanamori (G-K) rules 
|l7l give an indication of which exchange interactions 
should be positive (ferromagnetic) and which negative 
(antiferromagnetic) depending on the state of ionisation 
of the two ions, the occupied orbitals and the angle sub- 
tended at the bridging ion. They are valid only for in- 
sulating states and were worked out using perturbation 
theory to give a general guide to the interactions although 
deviations are known to occur. [T^ It is useful to compare 
our results for this specific material with the predictions 
of the G-K rules because the results may be used in fu- 
ture to assess the reliability of the rules. We note that 
in the case where we have ferromagnetism and ferromag- 
netically aligned orbits then our results predict a metallic 
ground state and so in these cases the rules are not ap- 
plicable. 

In LaMnOa the Mn ions are all in the same oxidation 
state and the Mn ions and the bridging oxygen lie along 
a straight line in the cubic unit cell; the bridging angle 
is 77. Thus the only variable that is relevant to the G-K 
rules is the orbital order. The rules state that if nearest 
neighbour sites are occupied by the same orbit the in- 
teraction is negative, antiferromagnetic. The size of the 
effect depends directly on the overlap of the orbits e.g. 
if there are two orbits — (which have large lobes 



in the x direction) separated by a lattice vector directed 
along X it would be larger than if, for example, the two 
orbits were — but still separated by a lattice vector 
directed along x. This would fit nicely with the result in 
Table ^ that the value of J2 (exchange up the z direc- 
tion) is negative for 00 — r'^ / — and also for 00 
— r^ /3z^ — r^ but significantly larger in the latter case 
where there are 3z^ — orbits arranged in columns up 
the z axis. The calculation of Ji is more complicated [13 
because the orbits are partially occupied but is ferromag- 
netic for 00 3a:^ — r^/3y^ — r^.(2Qj The overlaps would be 
smaller in the x-y plane for the case 00 3x^ — jiz^ — 
than for 00 Zx^ — jZy^ — so application of the G-K 
rules would predict a larger value of Ji (exchange in the 
x-y plane) in the former case in agreement with first prin- 
ciples results. The signs of Ji and J2 in Table ITU do agree 
with the G-K rules. There is one detail in which the first 
principles results do deviate from the G-K rules and that 
is in the case antiferromagnetic ordering with the ferro- 
magnetic orbital order F 3z^— r^. In this case since all the 
orbits are the same all the nearest neighbour interactions 
should be antiferromagnetic which would mean that the 
G-AF state should be more favourable than the A-AF 
state whereas the opposite order is seen in Table U and 
the value for Ji in Tabic ^1 should be negative for this 
orbit. The order is correct for the other ferromagnetic 
orbital order, F x^ — y^. Thus we see the predictions for 
the signs of Ji and J2 from the first principles calculation 
and the G-K rules agree in all cases (except the one men- 
tioned for F — orbit above) and the magnitudes 
agree. In one case we see that there is a disagreement on 
the ordering of unfavourable states. Model perturbation 
calculations of the exchange constants also disagree with 
the G-K rules: As mentioned earlier, Sikora and Oles 
9] have found that for the case oi 9 = 27r/3 the con- 
stants are small and both ferromagnetic, whereas G-K 
rules predict that Ji is strongly ferromagnetic while J2 
is antiferromagnetic. 

In the ferromagnetic case the total moment is 4/iB 
which is the value one expects from having four d elec- 
trons. This is the case also because the FM solution is 
either half-metallic (00 scenarios 3z^ — and x'^ — y^ of 
Table nj or insulating (00 scenarios 3x^ — r^/3y^ — r^, 
x^ - z^/y^ - z2 and Sx"^ - r'^/3z^ - of Tabled). The 
magnetic moment that is on the Mn ion can be less than 
this because of hybridiation . The magnetic moment on 
the Mn ion when one Cg orbit is localised is about 3.70/j,b 
in both FM and A-AFM solutions and of 3.60/i_B in the 
G-AFM case. Because of hybridisation with the oxygen 
part of the polarisation is sitting on the oxygen ion. 
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The system is insulating in both orbital ordering sce- 
narios independently of the magnetic ordering. Inspec- 
tion of the total density of states (DOS) in the lowest en- 
ergy — jy^ — ordering scenario presented in Figs. 
[TIA-AFM),|3FM) andElG-AGM) reveals the presence 
of a gap which is larger as more nearest neighbour spins 
become antiferromagnetic (See also Table . Its cal- 
culated value in the 3a;^ — jZy^ - orbital and AFM 
magnetic orderings is in very good agreement with the 
experimental optical gap '2l'| as can be seen in Table ITTTl 
The peak at about -0.75Ry in the total DOS corresponds 
to the localised "it^g and one Cg orbits. The latter are 
shown in Fig. 0] where we can see the following features 
in the majority spin channel: the peak at -0.75Ry rep- 
resenting the localised — states and the ix^ — 
states split into occupied states which hybridize strongly 
with the oxygen 2p states and unoccupied Zx^ — r^ states. 
One can also notice by looking at the minority eg states 
that both orbits are degenerate because these are not 
corrected for by the SIC and hence are solutions of the 
LSD potential which are orthogonal to the SIC states. 
In the LSD calculation the t^g and eg states lie near the 
Fermi level with the tig states somewhat more localised 
than the eg ones. However the LSD does not describe 
their localisation accurately. In the SIC they are pushed 
well below the valence band, composed mostly of Oxygen 
l-p states. It is however known that the position of the 
Sl-corrected levels does not correspond to what would be 
seen in experiment. Relaxation effects need to be consid- 
ered if one wanted to get spectra from SIC single particle 
energies Centred around -1.25Ry are the oxygen 2s 
and La 5p semi-core levels. 



TABLE III: Energy band gaps in eV. 



Configuration 


- rV3y^ - 


- z'ly' - 


Exp 


FM 


0.54 


0.27 




A-AFM 


1.09 


1.29 


1.1" 


G-AFM 


1.50 


1.56 





'Ref. ^ 



The total energy of the solution where only tig or- 
bits are localised and the eg electron is delocalised lies 
much higher than the most unfavourable orbital order- 
ing solution which confirms that there is strong tendency 
to the localisation of the eg electron in LaMnOa even 
in the cubic phase. The energy scale of the localisa- 
tion/delocalisation of the eg electron is indeed at least 
twice as big as the energy corresponding to ordering the 
orbits. This is qualitatively in agreement with the ex- 
perimental observation that even above the critical tem- 
perature of the orbital ordering local distortions remain. 
Local distortions are an indication that there is local- 
isation. Once these eg electrons are localised they in- 
duce local distortions through the interactions with the 
surrounding oxygens and these distortions order simul- 
taneously with the orbits when the temperature is low- 
ered. Although we can not with the current method sim- 



ulate real paramagnetism as being a collection of disor- 
dered local moments without long range ordering we can 
speculate however, since the orbital ordering is so strong 
and independent of the spin ordering, that orbital or- 
dering occurs in the paramagnetic state too. It is this 
orbital ordering which drives magnetic ordering and not 
the other way round. In a model calculation of param- 
agnetic LaMnOs and KCuFa based on an LDA-t-U elec- 
tronic structure Medvedeva et al. jl^ concluded that 
distortions were not needed to stabilise the orbitally or- 
dered phase in both compounds and that this ordering 
is of purely electronic origin. Their calculations for cu- 
bic LaMnOa have found that in the PM phase the orbits 
order but they are not pure local — and — y^- 
They found that the local Sz^ — has an occupancy 
of 0.81 and the local x'^ — y^ has an occupancy of 0.21. 
This is consistent with our present calculations in that 
the calculated ground state is nearly degenerate. In ear- 




Energy (Ry) 



FIG. 1: Total DOS in A-AFM magnetic and x'^ - z^ /y^ - z^ 
orbital orderings. 

lier LDA-I-U calculations "2^ on KCuFa it was found that 
within LDA there was no instability of the system against 
distortion while in LDA-I-U the energy has a minimum 
for a finite distortion of the lattice. It was concluded 
then that electron-phonon and exchange only are not 
enough to drive the collective distortion. A similar view 
was supported also by model calculations |3j US where 
both electron-electron and electron-lattice interaction are 
taken into account. In our present calculation the com- 
petition is rather in terms of localisation/delocalisation 
of the eg orbits by electronic interactions alone. And we 
found indeed that these are enough to first localise the 
orbits (larger energy scale) and then to order them in an 
anti- ferromagnetic way (smaller energy scale). Based on 
these results and those mentioned earlier we speculate 
that the distortions are a consequence of the displace- 
ment of oxygen ions to accommodate the electrostatic 
interactions resulting from the orbital ordering but these 
are crucial in selecting the ground state ordering out of 
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is still insulating even at higher pressure than the critical 
one at which the structural transition takes place. There 
was no indication of the magnetic state of the system but 
the experiments were carried out at room temperature 
which is well above the ordering temperature at least of 
the distorted LaMnOs. We found both FM and A-AFM 
solutions to be insulating in both Sx^ — r^/3y^ — and 
— /y^ — orbital ordered states. Whereas the sys- 
tem is metallic when only the t2g electrons are localised. 
The fact that the system was found to be insulating after 
suppression of the JT distortion is indicative of the pres- 
ence of orbital ordering with or without spin ordering. 
Use of local probe such as EXAFS or PDF (pair distribu- 
tion function) would be of great help though to settle the 
question of whether pressure really quenches distortions 
at the local level. 



FIG. 2: Total DOS in FM magnetic and ~ jy - - 
orbital orderings. 



the two nearly degenerate solutions we found for the cu- 
bic case. 




Energy (Ry) 

FIG. 3: Total DOS in G-AFM and 
orderings. 



/y^ — orbital 



Earlier SIC-LSD calculations by Tyer et al. ^ have 
described correctly the physics of the distorted LaMnOa. 
Then Banach and Temmerman |23| studied the cubic 
phase but using a unit cell of two formula units only. 
This limited the study to the first two rows and first four 
columns of Table U Hence they found that the lowest 
energy solution is the A-AFM with — orbital or- 
dering. Upon decreasing the lattice parameter they found 
a crossover to the FM with t2g orbitals Sl-corrected only 
which means suppression of orbital ordering. We recon- 
sidered this case below with our present bigger cell. 

Loa et al. studied structural and electronic proper- 
ties of LaMnOs under pressure and found that the system 



tn 
O 
Q 




-0.5 

Energy (Ry) 

FIG. 4: Mn eg-projected DOS in the ground state A-AFM 
in the — z'^ jy^ — 7? orbital ordering on the Mn site with 



the self interaction correction applied to the y 



orbital. 



Another way of suppressing the distortions is by in- 
creasing temperature as done by Sanchez et al. p9l | who 
studied the structural changes of LaMnOa with tempera- 
ture by using XANES and EXAFS measurements. Prob- 
ing the local environment of the Mn ions they found no 
abrupt change in the signal upon crossing the structural 
transition temperature Tjt- They described the struc- 
tural phase transition as ordering of the local distortions 
that are thermally disordered above Tjt resulting in a 
cubic lattice on average. This picture is quite different 
from the high pressure one although in both cases the 
distortions are apparently suppressed. In the high tem- 
perature regime orbital ordering can still be present but 
the long range ordering is suppressed by thermal fluctu- 
ations. Consistent with our calculation that the localisa- 
tion/delocalisation energy is of a larger scale than orbital 
ordering, i.e, the Cg electrons tend to localise strongly. 
As a consequence the lattice is distorted locally but since 
the energy scale of ordering the orbits/distortions is lower 
they are disordered by thermal fluctuations at high tem- 
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' A-AFM X -z /y -z 




-6 -4-2 2 

Lattice parameter(% deviation from a^^) 

FIG. 5: Total energies of the A-AFM - jy^ - orbital 
ordering and FM t^g solutions as functions of the deviation 
of the lattice parameter a from the experimental one ao = 
7.434a.u. 

state to the delocalised as a is decreased. 



perature. 

We have also investigated the dependence of the orbital 
ordering on the volume of LaMnOa. To do so we compare 
total energies for different lattice parameters relative to 
the experimental one. The latter is determined by re- 
quiring that it gives the correct experimental volume of 
the distorted system. We compared the energies of two 
scenarios: the ground state solution of the experimental 
volume (x^ — jy^ — orbital ordering and A-AFM 
spin ordering) and the FM solution with delocalised Cg 
orbits. The results are given in Fig. El One notices that 



the lattice parameter corresponding to the minimum is 
the same in both solutions and that it is slightly smaller 
than the parameter obtained from the experimental vol- 
ume of the distorted system. Upon decreasing the volume 
the two curves cross a at about -5% of the experimental 
lattice parameter. Below this value the electron be- 
comes delocalised and there is no longer orbital ordering. 
The system becomes metallic too as was signalled by the 
jump in the conductivity found by Loa et al. |28j . 



IV. 



CONCLUSIONS 



We have investigated orbital ordering in cubic LaMnOs 
using the SIC-LSD method which allows to study the 
localisation-delocalisation competition of correlated elec- 
trons. Although orbital ordering in LaMnOa has been 
ascribed to Jahn- Teller distortions of the MnOg octahe- 
dra we found that this ordering can happen from purely 
electronic effects by spontaneous breaking of the cubic 
symmetry. Once the orbital ordering sets in the electro- 
static interaction between the O ions and the electrons on 
the neighbouring Mn ions can be minimised by elongat- 
ing the bonds along the lobes of the occupied eg orbitals. 
It seems though that this coupling to the lattice is still 
needed to select the correct orbital ordering giving the 
observed distortions in the real LaMnOs system. There 
is therefore no need to assume an underlying A-AFM 
magnetic ordering to recover the orbital ordering. The 
latter is independent of the magnetic ordering and this 
is evidenced by the much higher ordering temperature of 
the orbits as compared to the spins. Although what we 
have found is that the lattice is important to determine 
the symmetry of the ground state orbital ordering. 



[1] Y. Murakami, J. P. Hill, D. Gibbs, M. Blume, L Koyama, 
M. Tanaka, H. Kawata, T. Arima, Y. Tokura, K. Hirota, 
and Y. Endoh, Phys. Rev. Lett. 81, 582(1998) 

[2] D. V. Efremov and D. I. Khomskii, Phys. Rev. B 72, 
012402(2005). 

[3] L V. Solovyev, N. Hamada, and K. Terakura, Phys. Rev. 

Lett. 76, 4825(1996). 
[4] W. E. Pickett and D. J. Singh, Phys. Rev. B 53, 

1146(1996). 

[5] S. Satpathy, Z. S. Popovic, and F. R. Vukaylovic, Phys. 

Rev. Lett. 76, 960(1996). 
[6] T. Mizokawa and A. Fujimori, Phys. Rev. B 51, 

12880(1995). 

[7] T. Mizokawa and A. Fujimori, Phys. Rev. B 54, 
5368(1996). 

[8] H. Zenia, G. A. Gehring, G. Banach, and W. M. Tern- 
merman Phys. Rev. B 71, 024416(2005). 

[9] O. Sikora and A. M. Oles, Acta Physica Polonica B 34, 
861(2002). 

[10] J. Deisenhofer, B. L Kochelaev, E. Shilova, A. M. Bal- 
bashov, A. Loidl, and H.-A Krug von Nidda, Phys. Rev. 



B 68, 214427(2003). 
[11] J. Rodriguez-Carvajal, M. Hennion, F. Moussa, A.H. 

Moudden, L. Pinsard, and A. Revcolevschi, Phys. Rev. 

B 57, R3189(1998). 
[12] J. P. Perdew and A. Zunger, Phys. Rev. B 23, 

5048(1981). 

[13] W. M. Temmerman, A. Svane, Z. Szotek, and H. Winter, 
in Electronic Density Functional Theory: Recent Progress 
and New Directions, edited by J. F. Dobson, G. Vignale, 
and M. P. Das, (Plenum, New York) 1998. 

[14] K. Hirota, N. Kaneko, A. Nishizawa, and Y. Endoh, J. 
Phys. Soc. Jpn. 65, 3736(1996). 

[15] F. Moussa, M. Hennion, J. Rodriguez-Carvajal, H. 
Moudden, L. Pinsard, and A. Revcolevschi, Phys. Rev. 
B 54, 15149(1996). 

[16] W. M. Temmerman, Z. Szotek, and H. Winter, Phys. 
Rev. B 47, 1184(1993). 

[17] J. B. Goodenough, Magnetism and Chemical Bond (In- 
terscience. New York, 1963). 

[18] D. Khomskii, Spin Electronics, edited by M. Ziese and 
M. J. Thornton (Springer, Berlin, 2001), p. 89. 



7 



[19] H Mcskinc, H Konig and S Satpathy, Phys Rev B 64, 
094433 (2001). 

[20] D.I. Khomskii and GA Sawatzky, Solid State Comm 102, 
87 (1997). 

[21] T. Arima, Y. Tokura, and J. B. Torrance, Phys. Rev. B 

48, 17006(1993). 
[22] J. E. Medvedeva, M. A. Korotin, V. I. Anisimov, and A. 

J. Freeman, Phys. Rev. B 65, 172413(2002). 
[23] A. I. Lichtenstein, V. I. Anisimov, and J. Zaanen, Phys. 

Rev. B 52, R5467(1995). 
[24] M. V. Mostovoy and D. I. Khomskii, Phys. Rev. Lett. 

92, 167201(2004). 
[25] S. Okamoto, S. Ishihara, S. Maekawa, Phys. Rev. B 65, 



144403(2002). 

[26] R. Tyor, W. M. Temmerman, Z. Szotk, G. Banach, A. 
Svanc, L. Petit, and G. A. Gehring, Europhys. Lett., 64, 
519(2004). 

[27] G. Banach and W. M. Temmerman, J. Phys.: Condens. 

Matter 16, 85633(2004). 
[28] I. Loa, P. Adler, A. Grzechnik, K. Syassen, U. Schwarz, 

M. Hanfland, G. Kh. Rozcnbcrg, P. Gorodctsky, and M. 

P. Pasternak, Phys. Rev. Lett. 87, 125501(2001). 
[29] M. C. Sanchez, G. Subias, J. Garcia, and J. Blasco, Phys. 

Rev. Lett. 90, 045503(2003). 



